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i 

Random  distributions  of  correlated  scatterers  averaged  over  orientation  are  considered, 
corresponding  to  isotropic  fluids  of  statistical  mechanics  particles  (with  volume  v,  number 
concentration  p,  and  volume  fraction  w  =  pv).  For  minimum  separation  of  centers  small 
compared  to  wavelength  and  acoustic  particle  parameters  close  to  the  embedding  medium’s, 
the  incoherent  differential  scattering  from  unit  volume  and  the  corresponding  attenuation 
coefficient  are  proportional  to  the  fluctuations  (variance)  in  number  concentration.  For 
arbitrary  convex  hard  particles  (e.g.,  ovals  or  simple  polyhedra,  repulsive  at  contact)  with 
shape  parameter  c>  3,  the  variance  is  expressed  in  terms  of  a  quotient  S(c;w )  of  polynomials  in 
w  that  has  a  maximum  S'r  (c)  at  w$  (c).  Spheres  (c  =  3)  were  considered  earlier.  For  c>  3, 
the  fluctuations  and  S*  and  w A  are  smaller  than  for  spheres;  for  c  <  3  (which  we  consider 
formally),  they  are  larger  than  for  spheres.  The  results  are  interpreted  by  comparing  leading 
terms  with  the  second  virial  coefficients  for  more  general  statistical  mechanics  models. 
Scattering  data  for  suspensions  of  diseoidal  red  blood  cells  versus  w  under  different  flow 
conditions  can  be  fitted  adequately  by  S(c;w)  for  different  values  of  c  <  3.  The  low  values  of  c 
suggest  weaker  repulsion  between  deformable  cells,  and  attractive  interparticle  forces  mediated 
by  flow  and  aggregative  trends.  ' 
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INTRODUCTION 

Earlier  papers1"4  analyzed  the  average  field  in  random 
distributions  of  correlated  scatterers  corresponding  to  fluids 
of  hard"  statistical  mechanics  particles  (with  volume  v, 
number  concentration p,  and  volume  fraction  w  —  pc).  The 
low  -frequency  scattering  and  the  corresponding  attenuation 
depend  on  the  fluctuations  (variance)  in  number  concentra¬ 
tion.  Using  the  scaled  particle*  approximate  equation  of 
state  E.  we  expressed  the  variance  in  terms  of  a  simple  quo¬ 
tient  S(u-)  of  polynomials  that  vanishes  at  the  extremes  ofw. 
For  spheres  (as  well  as  aligned  ellipsoids  ),S(w)  has  a  max- 
imumS  ~0.047atu’,  ~0. 129.  Now  we  consider  a  gener¬ 
alization  S(c.u')  for  isotropic  fluids  of  nonspherical  particles 
averaged  over  orientation,  with  c  as  a  parameter.  For  arbi¬ 
trary  convex  particles  (such  that  a  line  segment  connecting 
any  two  points  in  ris  wholly  within  tO.cis  determined  by  the 
volume,  surface,  and  average  over  angles  of  the  mean  of  the 
principal  radii  of  curvature-*  1(1  For  such  cases,  Gibbons" 
applied  scaled  particle  theory  in  terms  of  Isihara’s*  results 
for  two  different  convex  particles  to  obtain  Elc),  which  we 
use  to  construct  S(c;w). 

Convex  hard  particles  require  c>  3,  with  c  —  3  as  the 
special  case  of  spheres,  and  c>  3  as  a  nonsphericity  param¬ 
eter.  Fore  >  3,  we  show  that  SA  and  w  A  are  smaller  than  for 
spheres.  We  also  consider  3  >  c> 0  formally,  in  which  range 
they  are  larger:  0.047  5SA  50. 148  and  0.129  5  tuA  50.333. 

( The  Appendix  analyzes  E  and  5  for  c  <  0. )  The  behavior  of 
S(c;u> )  is  uniform  in  w  for  all  c> 0,  and  its  simplicity  suggests 


’’Work  supported  in  part  by  the  Office  of  Naval  Research  and  the  National 
Science  Foundation. 


inverting  data  showing  a  peak  at  0  <  u>  A  <  1  /3  to  obtain  an 
effective  c(wA  ).  For  c>  3,  we  expect  the  upper  bound  on  w 
to  be  less  than  for  spheres  (  s;0.63). 

For  a  hard  sphere,  the  cocentered  exclusion  volume 
( that  excludes  all  other  sphere  centers )  equals  8p.  The  aver¬ 
age  exclusion  volume  for  hard  ( repulsive  at  contact )  convex 
particles*  is  v,.(c)  —  2(c  -t-  l)u,  and  for  c>3  the  fluctu¬ 
ations  are  smaller  because  i\.  >  8 c  (less  elbow  room);  this 
also  holds  for  noncon vex  bodies "• 1 1  wit  h  an  effective  c  >  3 .  If 
we  use  the  form  v,.  (c)  for  c  <  3,  we  may  interpret  the  larger 
fluctuations  in  terms  of  an  effective  ve  <  8u  arising  from  ad¬ 
ditional  neighbors  at  small  separations;  comparison  of  lead¬ 
ing  terms  with  the  second  virial  coefficients14  for  nonhard 
particles  suggests  attractive  weakly  repulsive  models. 

Values  of  vL,  >  8tt  with  effective  c  >  3  do  not  imply  con¬ 
vexity;  see  Isihara’s  results"  for  nonconvex  particles  formed 
from  two  hard  spheres.  Although  the  work  by  Gibbons"  is 
based  on  convex  particles.  Rigby"  found  adequate  accord 
between  E(c)  for  the  nonconvex  tangent  dumbbell"  and 
Monte  Carlo  computations.  Similar  accord  is  shown 
between  S(c;w )  for  different  values  of  c  ^  3  and  reduced  ul¬ 
trasonic  backscattering  data15  versus  w  for  suspensions  of 
diseoidal  red  blood  cells  under  different  flow  conditions 
( stationary  or  stirred,  or  in  uniform  or  turbulent  flow ) .  The 
lower  effective  c  values  suggest  weaker  repulsion  between 
the  cells  (flexible  deformable  biconcave  discoids),  and  at¬ 
tractive  interparticle  forces  mediated  by  flow  and  aggrega¬ 
tive  trends;  values  nearer  3  may  arise  from  flow  alignment. 
Although  additional  factors  may  be  involved,  comparisons 
with  available  data 15  indicate  utility  of  S(c;w )  for  other  than 
hard  convex  particles. 
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For  minimum  separation  (b)  of  particle  centers  small 
compared  to  wavelength  (ir2/k)  and  acoustic  particle  pa¬ 
rameters  close  to  the  embedding  medium’s,  scattering 
aspects  are  particularly  simple.  The  fluctuation  function  S 
and  the  isolated  particle  scattering  (total  and  differential) 
ai.  J  absorption  cross  sections  suffice  for  the  coherent  attenu¬ 
ation  coefficient,  and  for  the  incoherent  scattering  from  unit 
volume  (except  for  translational  and  transducer  factors). 
Although  the  required  forms  can  be  obtained  by  simplifying 
earlier  results,  key  steps  of  direc'  derivations  are  included  to 
facilitate  applications  and  indicate  limitations.  We  also  con¬ 
sider  the  analogous  two-dimensional  problem  of  parallel  cyl¬ 
inders  and  obtain  the  corresponding  S(c;w)  by  using  a 
scaled  convex  disk1'’  E(c) ;  these  results  also  apply  for  mono- 
layers'*  of  bounded  particles.  The  one-dimensional  and  the 
lattice  gas  forms  of  5  are  included  for  comparison.  For  all 
cases,  the  function  E  corresponds  to  the  equation  of  state 
times  v  ( w  hich  leads  to  simpler  forms ) ,  but  we  use  no  special 
label  for  the  normalized  version. 


If  the  center  of  the  obstacle  is  at  r,  with  respect  to  the 
phaseorigin  (r  =  0),  then  we  replace  u ( r )  by  u{ r  —  r,  )e,kr'. 
For  large  r>r, ,  we  have  jr  —  r,  j  =/•  —  f*r, ,  and 

u(r  —  r,  )e,k'r'  =  u,A,  ~f(r,k)e'kik  n'r  (elkr/ r),  (7) 

where  A ,  =  <?,k  r’  is  the  phase  factor  introduced  by  the  inci¬ 
dent  wave. 

For  a  configuration  of  N  obstacles  in  a  volume  V,  the 
multiple  scattered  solution  has  the  form1’' 

q/  =  4S+  £  C/S(r  —  r, ).  (8) 

« -  -  i 

The  multiple  scattered  contribution  of  the  .sth  obstacle  (J, 
may  be  expressed  functionally  in  terms  of  the  single  scat¬ 
tered  amplitudes  of  all  obstacles  and  the  configurational 
variables  ( locations,  orientations,  etc. ) .  Multiplying  'F  by  an 
appropriate  probability  distribution  function  and  integrat¬ 
ing  over  all  variables,  we  write  the  average  wave  (the  coher¬ 
ent  field )  as 


I.  SCATTERING  ASPECTS 

We  consider  a  plane  wave  Ae  "" , 

A  =  e'11''.  k*r  =  Ark-r  =  kr  cos  6  =  kz,  ( 1 ) 

incident  on  an  acoustically  penetrable  obstacle  of  volume  t> 
with  center  at  r  =  0  (the  phase  origin  taken  as  the  center  of 
t he  smallest  sphere  enclosing  v ) .  The  corresponding  solution 
of  Helmholtz’s  equation  exterior  to  (has  the  form  A  =  A  u 
with  u  as  the  scattered  wave.  For  large  kr, 

u(r)  ~/(r,k)(e'*7r),  (2) 

where /(r,k)  is  the  conventionally  normalized  scattering 
amplitude  for  directions  of  incidence  k  and  observation  r. 
The  energy  transferred  via  interference  of  A  and  u  is  repre¬ 
sented  by 

Im  /fk,k)4,T/A  =  au  +  a,, 

(3) 

<?,  =  j1  /( r.k)  /  dll(r)  =  ^a-(f,k)rffl, 

where  aa  is  the  absorption  cross  section,  a,  is  the  scattering 
cross  section,  and  a{  r.k)  =  i  /( r.k )  i '  is  the  differential  scat¬ 
tering  cross  section. 

The  obstacle  is  specified  by  two  relative  acoustic  param¬ 
eters  C '  and  B ',  such  that  v':  =  C'/B'  with  r/'  as  the  relative 
index  of  refraction  in  v.  For  the  simplest  cases,  C'  is  the 
relative  compressibility  and  1/5 '  the  relative  mass  density; 
more  generally,  the  parameters  are  complex  with  Im  C'  >  0 
and  Im  6  '  <  0  to  account  for  absorption.  For  small  C'  ~  1 
and  5  —1,  and  largest  diameter  (2a)  small  compared  to 
wavelength,  we  use  Rayleigh's  form17 

/(r.k)  ~  [C'  —  1  —  (5  ’  —  I  )r*k]A  'v/\v  (4) 

for  arbitrary  shaped  v.  To  lowest  orders  in  k,  from  (3), 
oa  =  (Im  C  +  [Im  5'j  )kv,  (5) 

a(f,k)  =  |C’  -  1  -  (5'-  l)(rk)|2A  V/16jt, 

<r,  =  <|C'  -  l|2+  \B'  -  l|2/3)A  V/4jt,  (6) 

where  cr,  follows  by  integration,  essentially  as  in  Rayleigh’s 
developments. IK 


<VJ>>  =A  +  £<t/,>.  (9) 

The  difference  of  the  average  of  |'k]:  and  the  coherent  inten¬ 
sity  |('k>|2  is  the  incoherent  scattered  intensity 

/=  m-)  -  |<T>|: 

=  I<:^I:>  +z(z<W*>  -  (U, )£<[/?>).  (10) 

For  present  purposes,  we  need  consider  only  /  explicitly. 

We  replace  V,  in  ( 10)  by  the  single  scattered  wave  u,A ( 
with  implicit  averages  over  orientation,  etc.  Thus  the  single 
summation  reduces  to  fdr ,  p  with  p  =  N /V  as  the  average 
number  of  particles  in  unit  volume  for  the  homogeneous 
cases  of  interest.  The  double  summation  reduces  to 
p'Sdr,  fdr,  g(  |r,  —  r,|)  with  p~g  as  the  pair  distribution 
function  averaged  over  orientation.  In  terms  of  R  =  r,  —  r, 
for  the  separation  of  centers,  we  recast  the  integral  over  r,  as 
an  integral  over  R,  and  use  g(5)~  1  for  even  moderately 
large  R.  We  restrict  r,  to  a  relatively  small  central  region  Vc 
of  V,  a  region  containing  the  average  number  (n)  =  pVc  of 
irradiated  and  detected  particles,  and  work  with  the  farfield 
form  as  in  (7).  The  result 

/( f ,k )  =  p  y.  |  /( f ,k )  | 2  m  f  ,k )  /r , 


is  a  standard  form  for  x-ray  scattering  by  dense  gases  or 
liquids.2"  The  function  RTr.k)  is  the  statistical  mechanics 
structure  factor,  and  g  —  1  is  the  total  correlation  function. 
For  minimum  separation  of  centers  small  compared  to 
wavelength  (small  kb),  to  lowest  order  in  k. 


/(f,k)  =  4(r,k)f,/r, 
\(r,k)=pa(rM)W,  W=\+p[\g(R) 


(12) 
1  ]</R, 


where  W  is  the  low-frequency  limit  of  the  structure  factor. 

The  corresponding  attenuation  of  the  coherent  intensity 
|('I'(z)>|2  =  C(z)  for  single  path  propagation  of  (W)  in  V 
may  be  obtained  by  Rayleigh’s  procedure"'  for  the  ideal  gas 
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case.  Thus,  from  ( 12),  the  net  incoherent  scattering  from  the 
particles  in  unit  volume  is 


^(r.kjrfii  =pWos=  A,, 

(13) 

and  the  energy  they  absorb  equals 

r  =  ■ 

(14) 

The  sum 

T  =  +  A,  -p(otl  +  a,  W) 

(15) 

with  K  as  Boltzmann’s  constant,  T as  absolute  temperature, 
g  as  isothermal  compressibility,  and  p  as  fluid  pressure.  The 
final  equality  relates  W  to  a  derivative  of  the  equation  of  state 
p/KT  =  E  /v.  Using  the  scaled  particle"  forms  of  E  for  hard 
spheres,  circular  cylinders,  and  slabs  (E,  with  i  =  3,  2,  and 
1), 

„  w(  1  +  w  +  w2)  ..  w  „  w 

C  1  —  J  >  Cn  —  7  «  t*  I  —  t 

(1  -U’)  '  (1  -  ft’)'  1  -  w 

(22) 


provides  the  attenuation  coefficient  for  the  coherent  intensi¬ 
ty  in  the  form  dC(z)/dz  =  —  tC(z),  from  which 
C(z)  —  e  for  C(0)  =  1.  See  Ref.  2  for  full  development 
of  W. 

The  following  is  not  restricted  to  three-dimensional 
forms  of  W  and  the  single  particle  scattering  functions  in 
(6).  Analogs  for  two-dimensional  cases  of  normal  incidence 
on  parallel  cylinders  are  based  on 

rr(f.k)  =  C'  -  1  -  (B  -  l)r-k  :k  V78- 

,  ,  ,  .  06) 

rz,  =  ( i C  —  1  -  +  B‘  -  1  -/2)k  Y74. 

which  represent  cross  sections  per  unit  length  with  v  as  an 
area.  Similarly  for  one-dimensional  problems  of  normal  inci¬ 
dence  on  parallel  slabs. 

o(r.k)  =  C  -  1  —  ( B '  —  1  )rk  A  -’(••’/ 4. 

,  ,  ,  (H) 

<7,  -  (  C'  -  1  -  +  B'  -  I  -)A  r/2. 

which  represent  cross  sections  per  unit  area  with  r  as  slab 
thickness;  here  only  r  =  —  k  (backscattering)  and  r  =  k 
(forward  scattering)  apply.  The  corresponding  absorption 
cross  sections  (per  unit  length  or  unit  area)  are  given  by  (5) 
in  terms  of  the  present  r’s. 

II.  STATISTICAL  ASPECTS 

The  differential  incoherent  scattering  and  the  attenu¬ 
ation  via  scattering  are  governed  by  the  form 

A  =  p  Wo  =  So/v,  S  —  w  W.  w  — -  pi\  (18) 

where  o  is  either  the  differential  or  total  scattering  cross 
section,  and  w  is  the  volume  fraction  occupied  by  particles. 
The  packing  function  W.  determined  in  ( 12)  by  the  integral 
of  the  correlation  function  g  —  1,  is  proportional  to  the’" 
variance  v  (the  mean-square  deviation)  of  the  number  n  of 
particles  in  , 

v=  <(n  -  <«>)"'>  =  <«’)  -  (n):  =  (n)W.  (19) 

Thus 

S  =  u  W  =  \v/Vc  =  ({n:>  -  (n2)  )v/Vt .  (20) 

henceforth  the  fluctuation  function,  is  proportional  to  the 
mean-square  fluctuation  in  the  number  of  particles  around 
the  mean  (n)  —pVc. 

Under  appropriate  conditions,  a  distribution  of  particles 
may  be  regarded  as  a  large  scale  fluid.  The  variance  v,  and 
therefore  W,  can  then  be  related  to  statistical  thermodynam¬ 
ical  functions  by 

/  OTIC 


it  follows  from  w  —  pv  and  dE /dw  =  1/W  that 


Wt  = 


(I  -w)4 
(1  +  2 w)2  ’ 


W,  = 


(1  -  w)' 
1  +w 


W,  =  (  1  -w)2. 


(23) 


Form  E{  is  the  exact  Tonks’  equation  of  state,  and  Ex  was 
rederived  by:i  Wertheim  and  by  Thiele  from  the  Percus- 
Yevick  integral  equation  for g.  The  result  Wt  follows  directly 
by  integrating4  the  Zernicke-Prins  g,  and  both  W ,  and  Wt 
also  follow  from  the  Percus-Yevick  equation. 

The  packing  functions  W(w)  decrease  monotonically 
from  unity  to  zero  as  w  increases  from  zero  to  unity.  How¬ 
ever,  as  emphasized  before,1  the  fluctuation  functions 
S(w)  =  u’W(w)  vanish  at  the  extremes  of  icand  have  a  max¬ 
imum  S  (u’A)  at  w  =  w ,  corresponding  to  dS /aw 
=  S'  =  0  and  S  "  =  —  !5  " , .  Thus 


. .  w(  1  -  ic)4 
‘  (1  +  2u  ):  ’ 

(  ^4) 

tc,  -  [  (73) 1  -  -  71/12=0.129.  S.,  =0.0469; 

5  =  »’( 1  -  »’) ' 

1  +  (<  (25) 

=  (71-  -  2)/3  =  0.2 15.  S  =0.0856; 

2>,  =  w(  1  —  «’)’,  =1/3,  S.  =  4/27=0.148. 

(26) 


See  Fig.  1  for  plots  versus  u\  Although  we  may  consider  the 
full  range  0s.u’<l  for  some  purposes,  the  upper  bound  is 
physically  realizable  only  for  slabs.  For  spheres  and  cylin¬ 
ders  we  use  the  experimental  values  for  densest  random 
packing.  u\, ,  =  0.63  and  u\, ,  =  0. 84  ( corresponding  to  amor¬ 
phous  solids).  By  geometry,  the  results  for  spheres  also  ap¬ 
ply  for  aligned  ellipsoids,  and  the  results  for  circular  cylin¬ 
ders  also  apply  for  aligned  elliptic  cylinders.4 

In  addition  to  the  above  cases  /'  =  (3.2.1).  we  also 
showed'  that  the  symmetrical  function  S  =  u'(\  —  u') 
which  arose  in  an  earlier  development”  could  be  derived 
from  the  scaled  particle'’  E  for  a  random  lattice  gas: 

Eu  —  —  ln(l  —  w),  Wlt=l  —  w;  (27) 

5„  =  w(  1  —  it’),  wA  =1/2,  SA  =1/4.  (28) 

Shape  is  not  specified,  and  w  =  1  is  realizable  for  special 
processes  and  compliant  particles  that  pack  to  fill  all  space. 1 
Since  S„  is  not  a  special  case  of  the  closed  forms  derived  in 
Secs.  IV  and  V,  it  still  serves  for  less  correlated  contexts."” 
(See  Fig.  2.) 

The  fluctuation  function  determines  the  behavior  of  A 
vs  u\  but  the  total  attenuation  r  involves  an  additional  term  if 
absorption  is  present.  We  rewrite  (15)  as:> 
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FIG.  1  Plots  of  the  fluctuation  functions  S  of  (24) -(26)  versus  volume 
fraction  tv.  The  highest,  central,  and  lowest  curves  are  S,.  S:.  and  6\.  respec¬ 
tively.  (See  Fig.  21  of  Ref.  7  for  comparison  of  with  analogs  we  derived 
from  other  model' . )  In  order  to  show  complete  curve  shapes  ( here,  as  well 
as  in  all  subsequent  graphs  versus  ir).  the  full  ranged  U'-  l  is  displayed,  but 
the  physically  realizable  upper  bound  for  tr  is  smaller  in  general. 


i'T  —  o,w  -t-  cr.Sltr )  =  (<5ic  -}•  S(w)  ] <rs .  5  =  atl/as. 

(29) 

The  occurrence  of  a  maximum  r  atw  —  w  depends  on  the 
magnitude  of  8  and  the  requirements  w  =  ic  1  <  ic,, .  For 
small  8.  to  first  order  in  small  terms. 

rr  ~autc.  +o \S .  ,  w  ~  w  ■  +  0a/cT,  S  "■  !  >  w  ..  . 

(30) 


FIG.  2.  Plots  of  the  ideal  gas  function  5  =  w  (the  diagonal  line),  of  the 
lattice  gas  S„  ( the  simple  parabola )  of  ( 28 ) ,  and  of  the  set  shown  in  Fig  I  to 
delineate  the  relative  scales  of  the  fluctuation  effects  and  their  decrease  with 
increasing  exclusion  volume  (decreasing  elbow  room ). 


where,  eorresponding  to  i  =  (3,  2,  1,0),  the  second  deriva¬ 
tive  jS"  |  =  !S"(«’a  ) |  approximates  or  equals  (2.84,  2.21, 
2,  2).  Thus  absorption  shifts  the  maximum  of  r  to  larger  w. 
Because  8  is  small,  it  follows  that  for  sparse  concentrations 
the  attenuation  is  dominated  by  scattering  losses.  However, 
since  IV  decreases  with  increasing  w,  there  is  a  value  w  =  wx 
determined  W(u\  )  =  8  for  which  absorption  and  scattering 
losses  contribute  equally  to  the  attenuation  and  cr  —  2 au  u\ ; 
for  w>  wx  ,  absorption  dominates.  On  the  other  hand,  if  8  is 
not  small,  then  there  is  in  general  no  maximum  of  r  with 
variation  of  w.  If  8>  1,  then  absorption  dominates  the  at¬ 
tenuation  for  all  w. 


III.  SECOND  VIRIAL  COEFFICIENT 

To  facilitate  derivation  and  interpretation  of  the  gener¬ 
alizations  of  S  given  in  the  next  sections,  we  consider  the 
leading  terms  of  the  expansions  of  E  and  IV  in  powers  of  p 
(virial  series'4).  An  ideal  gas  distribution  corresponds  to 
E  =  pi'  and  W  =  1  (as  well  as  v  =  ( n )  and  S  =  w).  The  first 
departures  from  the  ideal  gas  values  are  proportional  to  the 
second  virial  coefficient14  /?,, 

E,/v,  =  p  +  pB>  +  •  •  •  .  W,  =  1  -  p2B:  +  ■■■  , 


with  i\.  =  2  ‘  v,  and  r,  =  (t-4uV3.  -  a'.  2a.  r).  The  second 
term  of  W,  for  /  =  (3,  2.  1 )  also  follows  directly  from  (12) 
by  using  the  hole  approximation  for  the  radial  distribution 
function:  g(R )  =  0  for  R  <  b,  and  g(.R)  =  1  for  R  b.  with 
6  -  2o  as  the  minimum  separation  of  particle  centers. 

The  size  of  the  normalized  exclusion  volume  v,  /i\ 
~  21  —  (8,4,2,  1 )  determines  the  departure  from  the  ideal 
gas  results.  Although  ( 3 1 )  is  restricted  in  general  to  small  w. 
comparison  with  5,  shows  that  as  i\./v,  increases,  the  fluctu¬ 
ations  decrease  for  all  w  (less  available  space);  the  peak  S r, 
decreases  and  its  location  u>A  shifts  to  smaller  values 
(sparser  concentrations).  As  the  ratio  i\./v,  doubles  in  size, 
the  ratios  of  ajacent  pairs  of  5,  A  values  (0.55,0.57, 0.59)  and 
of  adjacent  pairs  of  «’,A  values  (0.60,  0.645,  0.67)  are  of  the 
order  of  half.  Thus  the  second  virial  coefficient  B:  not  only 
indicates  nonideal  gas  behavior  at  small  w,  but  also  supports 
major  trends  of  the  maximum  fluctuations  SA  and  of  their 
locations  wA  .  The  generalizations  of  5  derived  in  subse¬ 
quent  sections  exhibit  similar  relations  with  B:. 

Tim  exact  second  virial  coefficient  B2  for  an  isotropic  gas 
of  a  ‘  cl  i. my  convex  hard  particles  was  derived  by  Isihara,K 
who  ,t  nined  the  exclusion  volume  of  a  fixed  particle  by 
moving  another  in  contact  around  it  at  fixed  relative  orienta¬ 
tion,  and  then  averaging  over  orientation.  The  resulting 
aver|ge  exclusion  region  i\.  =  2 B2  equals 


tV  =  2 (i>  +  sr)  —  2v(  1  -I-  c), 

sr  _  f  d(l 

c  =  — ,  r-  (r,  +  r.) - , 

V  J  877 


(32) 


with  v  and  s  as  the  volume  and  surface  area  of  the  particle, 
and  7  as  the  average  over  all  angles  of  the  mean  of  the  parti¬ 
cle's  principal  radii  of  curvature  ( r ,  +  r:)/2.  For  spheres. 
c  =  3  and  =  8c  as  before;  for  all  other  bounded  convex 
particles,  c>  3andc(,  >  8c.  Isihara’s  result  for  c,  (and  for  the 
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average  exclusion  volume  for  two  different  convex  particles) 
was  obtained  independently  by  Hadwiger,1'  and  their  publi¬ 
cations  and  those  of  Kihara"’  provide  additional  relations 
and  various  illustrations. 

The  nonsphericity  parameter  c  >  3  in  ( 32 )  does  not  rep¬ 
resent  a  unique  shape.  For  example,  for  spheroids  (both  ob¬ 
late  and  prolate),  Isihara*  obtained  the  single  form 


w  ith  a  and  a'  as  the  largest  and  smallest  semidiameters.  Thus 
for  the  same  value  of  a'/a.  both  the  oblate  and  the  prolate 
shapes  give  the  same  value  for  c  (corresponding  essentially 
to  an  interchange  of  the  forms  s  and  7.  and  using  the  appro¬ 
priate  t'Y).  For  near  spheres,  a~a',  and* 1,1 

f;3  +  4(e4  +  0/15.  (34) 

For  near  disks  and  near  needles,  a ’pa',  and 

c  ~  3/4  +  ir3a/a'%  (35) 

is  w  ithin  4T  of  (  33 )  for  a/a'  $  3.  In  this  range,  c  increases 
linearly  with  a/a',  in  accord  with  Isihara's  almost  linear 
curve  based  on  (33).  From  the  complete  form  (33),  we  have 

a/a'  =  (6. 5. 4.3, 2. 1. 5), 

(36) 

r;  ( 7.65,6.55,5.48.4.45.3.54.3.18), 

w  here  c  decreases  to  the  sphere  value  3  for  a/a'  =  1. 

For  nonsmoo’h  convex  particles,  front  HadwigerY'  val¬ 
ues  (ofc.  v.  and  7 )  for  the  regular  polyhedra  ofy  faces,  we 
construct 

j  -  (4.6.8, 12.20).  r;  (6.7,4.5,4.32.3.56,3  45).  (37) 

The  largest  c  corresponds  to  tetrahedra;  the  progression 
through  cubes,  octahedra.  dodecahedra,  ends  with  icosahe- 
dra  with  <  only  1597  larger  than  for  spheres. 

It  is  clear  that  the  value  of  the  nonsphericity  parameter 
for  any  particular  polyhedron  can  be  matched  by  that  for  a 
spheroid  pair.  Thus  the  tetrahedron  value  cs6.7  corre¬ 
sponds  to  spheroids  with  a/a  zz  5.14:  the  markedly  jagged 
polyhedron  is  equivalent  to  a  markedly  flattened  or  elongat¬ 
ed  spheroid,  and  all  three  shapes  yield  v,./vzz.  15.4,  about 
tw  ice  the  value  for  spheres.  (The  closed  form  S  of  the  next 
section  shows  that  the  corresponding  values  of  S A  and  to  A 
are  of  the  order  of  half  the  sphere  values,  as  may  be  inferred 
from  the  discussion  following  (31 ).]  This  aspect  is  stressed 
to  indicate  the  limitations  of  c  as  a  shape  parameter;  we  may 
invert  data  in  order  to  obtain  c  from  measurements  of  B.  for 
particles  having  the  same  volume  c.  but  more  than  c  is  re¬ 
quired  to  determine  particle  shape.  For  that  matter,  an  effec¬ 
tive  c  >  3  isolated  by  measurements  need  not  correspond  to  a 
convex  particle  (or  to  a  hard  particle). 

Existing  results  for  B2  for  other  than  hard  convex  parti¬ 
cles,1’  14  can  be  written  formally  as 

B:  =  v,./2  =  u(c  +  1 )  (38) 

to  obtain  associated  values  of  ve  andc.  Although  c  =sr/v> 3 
is  an  exact  relation  for  convex  hard  particles,  an  effective 
c>3  in  (38)  does  not  imply  convexity.  In  particular,  Isi- 


hara 12  obtained  the  exact  B2  for  the  infinite  set  of  nonconvex 
hard  particles  (dumbbells)  formed  from  two  identical 
spheres  for  all  values  of  the  center  separation  b>  0.  (Only 
the  degenerate  case  b  =  0  for  a  single  sphere  corresponds  to  a 
convex  particle,  c  =  3.)  From  his  results,12  we  write 


b /a  —  (0,1,2,3,4,00  ),  c  =  (3 
b /a>4 ,  c  =  7  —  16(a/fi)2/5 


£■=  (3,3.329,4.444,6.031,6.8,7), 

(3Q) 


The  parameter  c  increases  smoothly  from  the  lower  bound  3 
to  the  upper  bound  7  (for  which  each  dumbbell  reduces  to 
two  independent  spheres). 

We  can  obtain  effective  values  c<3  (as  well  as  c>3) 
from  (38)  by  comparison  with  forms  of  B:  for  more  general 
models  than  hard  particles.  Although  hard  particles  exert  no 
forces  on  each  other  except  repulsion  at  contact  ( pair  poten¬ 
tial  infinite  for  R%b,  and  zero  for  R>b),  other  statistical 
mechanics  models14  include  attractive  forces  for  small  sepa¬ 
ration  (square  well  potential,  negative  for  b  <,Rc.  Ab).  as 
well  as  longer  range  weaker  repulsion  and  smoother  attrac¬ 
tion  ( Lennard-Jones  potential,  etc.).  Analogs  for  the  pres¬ 
ent  development  can  be  based  on  the  integral  in  (12)  with 
g(R  )  independent  ofp.  In  particular,  the  square  well  poten¬ 
tial  corresponds  to  supplementing  the  hole  approximation, 
discussed  after  (31).  with  g{R)  -  1  =  T  for  b  <  R  I b 
( representing  additional  neighbors  at  small  separation )  to 
construct 

IB.  =  2'r,  [1  -  (.4  '  —  1 ) r,  ] .  /=  (3.2.1).  (40) 

For  T,  --  0.  we  have  2'r,  as  before:  we  may  also  obtain  the 
lattice  gas  result  by  using  .4  =  2  and  T,  -  1/2'.  For  spheres. 
from  (40)  and  (38). 

2 B,  =  8 r(  1  -  (.-I1  -  1  )T]  -  2r(  1  +  c). 

<•  =  3  - 4(,4  1  -  nr.  (4I) 

so  that,  e  g.,  if  .4  =  2,  then  c—  (2.1.0)  corresponds  to 
28T  —  (  1.2.3 ).  etc.  The  size  of  T  can  be  related  to  an  effec¬ 
tive  interparticle  attractive  potential,  and  the  size  of  .4  to  the 
range  over  which  it  acts.  ( Changing  the  sign  of  V  yields  c  >  3 
corresponding  to  fewer  neighbors  and  repulsion  for 
b  <  R-  ,4b. )  Equation  (41 )  serves  for  interpretation,  but  no 
closed  form  of  5  is  available  for  this  model. 

For  two-dimensional  hard  convex  particles. we  write 

r  =  2c  f  xr  —  r( 2  +  c).  c  =  sr/v  -  s:/2~v. 

r  (42) 

v  =  I  r,  dO  --  2rrr. 

with  v  as  (he  cross-sectional  area  (or  volume  per  unit 
length),  v  as  the  length  of  the  perimeter,  and  7  as  the  average 
of  the  radius  of  curvature  ( r, ) .  For  circles.  c  =  2  and 
v,.  =  4c;  for  all  other  two-dimensional  convex  hard  figures, 
c  >2  and  c,.  >  4c.  As  before,  a  particular  value  of  c  does  not 
represent  a  unique  shape. 

For  ellipses. 


c  =  ( Sa/o 


/•ir/2 

V2)^2,  /  — 

Jo 


(1  —  e2  sin2  0) 1 


with  /  as  the  complete  elliptic  integral  of  the  second  kind. 
For  near  circles. 


:  ( 2a /a  )  ( 


eJ/32). 
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(46) 


(48) 


and,  for  near  strips, 

c  =  (8u/a';r){l  +  [ln(4<j/a')  —  1  ]  (a'/a)2},  af>a'.  (45) 
From  (43),  the  values  for  a  set  of  ellipses 
a/a'  =  (6, 5,4, 3,2, 1.5), 
c  =  ( 5.24,4.47,3.73,3.07,2.38,2. 13), 

decrease  to  the  circle  value  c  —  2  for  a/a'  =  1. 

For  nonsmooth  shapes,  we  consider  y-sided  regular 
polygons  in  terms  of 

c  =  (2y/7r)tan  (ir/j),  (47) 

and  obtain 

j  =  (3.4,5,6,8,12), 
c  =  (3.31,2.55,2.31,2.2,2. 1 1,2.05). 

The  triangle  value  c  =  3.3  also  arises  for  the  ellipse  with  a/ 
a'  —  3.45;  the  corresponding  vj i’  =  5.3  is  about  33%  larger 
than  for  circles. 

From  (38)  and  (40).  the  two-dimensional  square  well 
analog  of  (41 )  is 

25,  =4c[  1  -  (A  -  -  nr]  =  i  (2  +  c). 
c  =  2  -  4(a  2  -  nr 

so  that,  e  g.,  if  .4  =  2,  then  c*  =  (1,0)  corresponds  to 
1 2 r  =  ( 1,2),  etc. 

IV.  GENERALIZED  THREE-DIMENSIONAL  S 

Gibbons"  applied  scaled  particle  theory''  in  terms  of 
Isihara's  result'1  for  the  average  exclusion  volume  for  two 
different  particles  to  obtain 

E  =  tf|  1  4-  w(c  —  2)  -c  w:(  1  -  c  +  c2/ 3)  ]/( 1  —  w) 

(50) 

From  dE  /du'  =  E'  =  \/W,  we  construct 

id  1  -  if)4 


(49) 


S(c,ic)  =  u’W'(c;«')  = 


(51) 


[I  4  (C  —  1  )w 

The  present  E  and  5  correspond  to  convex  hard  particles  if 
c  ,-3  as  in  ( 32).  and  reduce  to  the  sphere  results  £\  and  5,  for 
e  -  3.  To  0{u,:)  we  obtain  (31 )  in  terms  of  (32).  At  u> 
=  1/2. 

S(c;  1/2 )  =  l/8(c  +  1):  =  (ivV..  );/2  (52) 

exhibits  the  direct  dependence  on  the  normalized  exclusion 
volume. 

The  fluctuation  function  S(c.u')  is  stationary  if 
3(r  -  I)tr’ +  (4+  f)«>-  I  =0  (53) 

corresponding  to  dS /dw  =  S'  =  0,  and  has  a  maximum 
Sic.w  ■  )  at 

( 4  +  20c  +  c: ) 


it’.  ( c )  = 


(4  +  c) 


6(c  —  1 ) 
2 


(54) 


4  +  c  +  (4  +  20c  4-  c:) 1  2 
From  ( 53 )  we  also  have 

c  =  1  +  ( 1  -  5 ic .  )/u\  ( 1  +  iwA  )  (55) 

and  eliminate  c  from  S(c;itv  )  =  SA  to  express  the  maxi¬ 
mum 

J  Acoust  Soc  Am  ,  Vol.  81 .  No  5.  May  1 987 


FIG.  3  Plots  of  the  generalized  three-dimensional  fluctuation  function 
5(c;u')  of  (51 )  with  the  parameter c  ranging  from  0  ( the  highest  curve )  to  8 
(the  lowest)  in  steps  of  1.  The  darker  curve  (r  -  3)  is  the  same  as  for 
hard  spheres  (or  aligned  ellipsoids);  the  lower  curves  ( c  >  3 )  apply  for 
( 59).  and  the  higher  curves  (e<  3)  for  (60)  with  the  highest  (c  0)as5, 
The  maxima  (the  peaks  S  )  and  their  locations  ( ir  )  on  the  w  scale  are 
plotted  in  Figs.  4-6. 


5-  =  ic.  (1 


)-(l  3ir  •  )-/4 


(56) 


solely  in  terms  of  its  location  on  the  volume  fraction  axis.  We 
include 


5  7  = 


d:S(c,w )  I 
du'  I  „  „ 

(!-«'•  )'[5  +  (c-  l)(l  +  6ic  .  )] 
ii  +  (c-  nicA  r 
(1  +  3u'a  ):( 1  +  6 ,  —  15ic:.  ) 

8ir  • 


(57) 


tofacilitate  applications  of  (30).  See  Figs.  3-6  for  plots  of  all 
key  aspects  of  (51  )-(56). 

We  can  apply  the  above  to  all  convex  hard  particles, 
including  the  spheroids  and  polyhedra  considered  in  (36) 
and  (37),  subject  only  to  the  restrictions  of  the  scaled  parti¬ 
cle  development.  For  example,  for  spheroids  with  a/a'  =  4. 


c=:5.478,  u-A  -0.0932,  SA  =0.0314, 


S":  =  -  3.14. 

(58) 


The  maximum  for  spheres  as  in  (24)  is  about  50%  larger 
and  occurs  at  a  packing  fraction  about  38%  larger;  the  pres¬ 
ent  —  5"  is  about  10%  larger  than  the  sphere  value  2.84. 
As  indicated  in  Sec.  Ill,  the  differences  arise  from  the  larger 
size  of  the  exclusion  volume  for  spheroids  averaged  over  ori¬ 
entation:  The  present  cr  /v ~  13  is  about  62%  larger  than  8 
for  spheres,  and  there  is  less  room  for  fluctuations.  As  men¬ 
tioned  after  (26),  the  sphere  form  S,  also  holds  for  aligned 
ellipsoids;  thus  the  above  also  serves  for  comparisons  of 
aligned  and  randomly  oriented  spheroids. 
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FIG  4.  Peak  5.  vs  c,  based  on  (51 )  in  terms  of  (54),  corresponding  to 
Fig  3. 


in 


Fig.  3. 


Although  S(c:iv)  has  the  same  value  for  both  the  oblate 
(4,  4.  1  )a'  and  prolate  (4,  1.  1  )a‘  cases  for  given  u\  the 
prolate  corresponds  to  four  times  as  many  particles  in  unit 
volume;  for  the  same  packing  fraction,  both  (n)  —  pVc.  (the 
average  number  in  the  central  region  \\  )  and  the  variance 
i  =  (/i)  IV  =  Sy  /r  are  four  times  larger  for  the  prolate 
case.  However,  the  normalized  scattering  cross  sections 
A  =  Sa/v>xSt  (and  particle  volumes)  are  four  times  larger 
for  the  oblate  case.  For  special  processes,  the  present  isotrop¬ 
ic  results  may  be  restricted  to  moderate  values  of  w,  with 
increasing  packing  the  particles  may  impede  each  others  ro¬ 
tational  freedom,  and  transitions  to  aligned  distributions 
may  arise.  For  such  processes  involving  hard  ellipsoids,  al- 


O  . 
o  1 

i  i  I - 1 - 1 - 1 

3  2  «  6  8  10 


FIG  5  Peak  location  w  vs  r.  based  on  (54),  corresponding  lo  Fig  .V 
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though  we  expect  S{c:w)  and  S ,(«■)  to  suffice  for  small  and 
large  w.  respectively,  the  transition  region  need  not  be  uni¬ 
form  and  could  entail  a  phase  change.  [For  spheroids  aver¬ 
aged  over  orientation  within  acoustically  transparent 
spheres'  of  radius  4  a'  and  volume  r„  we  have  A  =  Sx(u\)cr/ 
r,  or  A>7t’„  so  that  the  scattering  is  sixteen  times  larger  for 
the  oblate  than  for  the  prolate  case.  The  present  scattering 
maxima  A  .  ,.  and  A  •  ,,  are  about  2.68  and  10.7  times  larger 
than  the  sphere  constrained  analogs,  and  occur  at  w  about 
72 %  of  the  sphere  value.] 

Although  the  derivation  of  E  and  5  as  in  ( 50)  and  (  51 ) 
is  based  on  convex  hard  particles,  the  simplicity  of  the  results 
suggests  heuristic  applications  to  other  cases  for  which  no 
closed  forms  in  tr  arc  available.  Rigby 11  used  E  of  (50)  and 
Isihara's1'  value  for  the  nonconvex  contact  dumbbell 
(b  =  2a)  with  equivalent  c as  in  (39),  and  reported  satisfac¬ 
tory  accord  with  Monte  Carlo  computations.  Similar  accord 
is  shown  between  A  of  (51)  for  different  values  of  c  $  3  and 
backscattering  data1'  for  suspensions  of  particles  under  dif¬ 
ferent  flow  conditions.  We  therefore  consider  S(c.'r)  for¬ 
mally  for  allr>0:  In  this  range,  Sis  nonsingular  for  all  w,  and 
has  a  maximum  SA  <4/27  corresponding  to  w A  <  1/3  as  in 
( 54 )- ( 57 )  The  bounds  are  attained  for  c  =  0,  which  repro¬ 
duces  S,  of  (26).  (The  range  c  <0  is  analyzed  in  the  Appen¬ 
dix.) 

To  display  the  dependence  of  u<A  and  SA  on  c,  we  in¬ 
clude  two  sets  of  examples,  one  for  c>3,  and  the  other  for 
3>c>0.  Thus,  covering  the  illustrations  in  (36)  and  (37), 
and  including  (39)  formally, 

c=  (8,7,6,5.4,3), 

I00u\  cr (7.38,8.03,8.83,9.82,1 1.1.12.9).  (59) 

100SA  =  (2.36,2.62,2.94,3.35.3.9,4.69). 

where  the  final  entries  correspond  to  S,.  The  second  set  is 
formal,  with  interpretations  suggested  by  (41 )  ; 
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(60) 


c  =  (2.5,2,1.5,1,0.5,0), 

lOOte A  -  ( 14, 1 5.5, 17.4,20,24.2,3.33 ), 

10OSA  -(5.23,5.92,6.9,8.2,10.3,14.8), 


where  the  final  entries  correspond  to  S’,.  The  essential 
aspects  of  the  two  sets,  and  of  the  general  form  S(c;w)  of 
(51 ),  change  uniformly  with  c  in  the  range  c> 0. 

For  very  large  c,  we  have  w  A  —  1  /  (c  +  7 ) .  The  general¬ 
ization 

teA  ~  [4  4-  c  +  3(c  —  1  )/(4  4-  c)  ]  1  (61) 


is  within  2 %  of  (54)  for  c>  1.  At  c  =  1,  we  obtain  a  simple 
polynomial  for  5, 

c  —  1 ,  if,  =0.2,  5,  =0.082,  S  =  u'(l -tf)J.  (62) 


For  c  near  0, 

if,  -(1  -  c  4-  2c:)/3,  (63) 

and  5  reduces  to  the  simple  polynomial  5,  of  (26)  at  c  =  0. 


V.  GENERALIZED  TWO-DIMENSIONAL  S 

Boublik"’  applied  scaled  particle  theory"  in  terms  of  re¬ 
sults  for  two  different  hard  convex  disks  to  obtain 


r  !f[l  -ruffe-  2)/2] 

E’ - HTTP - ■ 

Front  E  '  =  1/W,  we  construct 


(64) 


S(c:w )  =  uH(c;w)  =  — — - — —  .  (65) 

1  T  (c  —  I  )ic 

The  present  E  and  S  correspond  to  convex  hard  figures  if 
c  2  as  in  (42 ),  and  reduce  to  the  circle  results  E  and  5,  for 
c  2.  To  O(w')  we  obtain  (31)  in  terms  of  (42).  At  w 
1/3. 


S(r.  1/3 )  =  8/27(2  4-c)  =  8f/f.  27  (60) 

exhibits  the  direct  dependence  on  the  exclusion  region. 

The  fluctuation  function  is  stationary  ( S'  -  0)  if 

3(c  1  )u"  -  4tf  -1=0  (67) 

and  has  a  maximum  S  at 


3tf  -  I;  2  +  ( i  -t-  3e ) ' 


Using 

e  =  I  4  ( 1  -  4tf  .  )/3tf  (69) 

we  eliminate  c  from  5  ■  =  S(c:tc  .  )  to  obtain 

5.  =  3tc:.  ( 1  -  if  ,  )’.  (70) 

The  corresponding  second  derivative 
S  "  =  --  ( I  -  if  .  ):[4  +  6(c-  I  )w  .  )/ 

[  1  4  (f  -  1  )u\,  |  =  -  1 8tf  .  (  1  -  2tf  A  )  (71 ) 


is  required  for  ( 30).  See  Figs.  7- 10  for  plots  of  all  key  aspects 
of  (65) -(70). 

We  can  apply  the  above  to  any  of  the  cases  in  (46)  and 
( 48 ),  e  g.,  for  elliptic  cylinders  with  a/a'  =  4, 

c  —  3.728,  ic.  -0.182,  S,  -0.0666,  S"  s  -  2.084. 

(72) 


l r> 


FIG.  7.  The  generalized  two-dimensional  fluctuation  function  5(c;«  )  of 
(65)  with  the  parameter  c  ranging  from  0  (the  highest  curve)  to  6  (the 
lowest)  in  steps  of  1.  The  darker  curve  <c  -  2)  is  the  same  as  .V.  for  hard 
circular  cylinders  (or  aligned  elliptic  cylinders);  the  lower  curves  (c  .  2) 
apply  for  (73).  and  the  higher  curves  (c*-  2)  for  (74)  with  the  highest 
<c  0)  as  6‘,. 


The  maximum  for  circles  (or  aligned  ellipses4)  as  in  (25)  is 
about  29<T  larger  and  occurs  at  a  packing  fraction  about 
18 cf(  larger.  |For  the  constrained  average  over  orientation 
within  acoustically  transparent  circles*  of  radius  4a'  and  vol¬ 
ume  e„  we  use  S:  of  (25)  and  A  =  S:a/r:  a  S,(u’:)r:/r: 
with  r ,  =  4f.  The  present  scattering  maximum  A  • 
=  S,  a/v  is  about  3.1  times  larger  than  for  the  circle  con¬ 
strained  maximum,  and  occurs  at  if  about  85°?-  of  the  cir¬ 
cle  value.  ] 


FIG  X.  Peak  5  vs  r.  based  on  (65)  in  terms  of  (68).  corresponding  to 
Fig.  7. 
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4  6  8  10 


FIG.  Peak  location  w  vs  c.  based  on  (68).  corresponding  to  Fig.  7. 

We  proceed  as  for  (51 )  and  consider  (65)  for  all  c  -0. 
Thus,  covering  the  illustrations  in  (46)  and  (48).  and  more 
generally  on  a  formal  basis  for  c  >  2, 

c  =  (7. 6,5, 4,3. 2). 

lOOtr.  -  ( 1 .49,  15.7,  16.7.  17.8,  16.4,21.5).  (73) 

1005.  -  (4.85.  5.27.  5.76,  6.44,  7.32,  8.56). 

where  the  final  entries  correspond  to  5:.  For  2  >c>0,  for¬ 
mally.  with  interpretations  suggested  by  (46), 

cm  (1.5,  1,0.75.0.5,0.25.0). 

100tc ,  -  (23,  25,  26.3,  27.6,  30.1.  33.3),  (74) 

1005.  s(9.42,  10.5,  11.3,  12.2,  13.3,  14.8). 


U“> 


C.C0  C.C5  0.10  0.15  0.2C  0.25  C.30  0.35 


W  RT  DERK 

HG  10  Peak  S  \s  u  tivithr  implicit ),  based  on  ( 70),  corresponding  to 
Fig.  7. 


where  the  final  entries  correspond  to  5,,  as  before  for  (60) 
The  essential  aspects  of  the  two  sets,  and  of5(c;n>)  of  (65), 
change  uniformly  withe  in  the  range  c>0.  (See  Appendix  for 

c<0.) 

For  very  large  c,  we  have  wA  —  l/v3c.  More  generally, 
we  use  the  complete  form  in  ( 68 ) ,  but  approximations  for  c 
near  1  and  0  are  of  interest.  Thus,  for  c  —  1, 

u>A  =(4  +  3(c-  1)/41  '  •  (75) 

At  c  =  1,  we  obtain  the  analog  of  (62), 
c=l.  w.  =0.25.  5,-0.105.  5  -  U’(  1  —«')'.  (76) 
For  c  -0, 

wA  =(1  -c/2  +  5c78)/3  (77) 

and  5  reduces  to  5,  at  c  =  0. 

APPENDIX 

For  c  <  0,  in  terms  of  t  =  c  —  1  and  d  =  —  t  >  0  (for 
brevity).  5  is  singular  at  u\,  =  \/d.  In  three  dimensions,  E  is 
nondecreasing  and  5  is  non-negative  for  all  w.  and  we  consid¬ 
er  the  full  range  for  insight.  In  two  dimensions,  E  decreases 
and  5 becomes  negative  for  w  >  u\,.  and  we  restrict  consider¬ 
ation  to  w  <  u\r 

We  rewrite  the  three-dimensional  equation  of  state  E  of 
( 50 )  as 

E  =  tc|  1  +  tc(7  —  1 )  +  w:(t 2  -  /  +  I  )/3]/(  1  —  tc) 1 , 

(Al) 

which  corresponds  to  hard  convex  particles  for  /> 2.  Here. 

I  =  2  gives  Ex  for  spheres,  and  t  =  —  1  gives  Et  for  the 
bound  c  =  0  considered  in  the  text.  The  first  two  derivatives 
with  respect  to  «•  are 

£  ’  —  t  '  4-  !lv)~  _  1 

“  (1  -«’)J  ~  w' 

( A2 ) 

2(2  -t - 1  +  tu<){  1  +  tw)  W' 

(1  -  te)'  “  ~lv2' 

If  t>  —  1,  then  E  increases  monotonically  to  oc  as  w  in¬ 
creases  tol.If  -r  =  d>\ then  at  u\,  —  —  I  //  =  \/d  there 
is  a  point  of  inflection  (£'  =  £"=  0)  for  which  E  =  w/ 
3(  1  —  w)  =  1/3  (d  —  1 );  however  (for  any  real  value  of  I), 
E  cannot  vanish  for  w  >  0,  and  is  nondecreasing  for  increas¬ 
ing  w.  The  point  of  inflection  suggests  a  phase  transition,  and 
the  corresponding  singularity  of  fV  and  5  =  w  W  at  wd  sug¬ 
gests  the  critical  region  (infinite  compressibility  C ). 

The  fluctuation  function 

5  =  icW  =  w(  1  —  w)4/(  1  +  tw)2  (A3) 

is  stationary  (5'  =  0)  at 

w  ,  =  {  -  (5  +  t)  ±  [(5 +  /):+  12/  ;,  :}/6t .  ( A4) 
Both  roots  are  of  interest  for  t  =  —  d,  as  are  both  values  of 
5 "  =  —  ( 1  —  w)  '(5  +  t  +  btw)/(  1  +  tw) ',  w  =  w  ,  . 

( A5) 

For  all  —  t  =  d>  1,  we  see  that  5  is  singular  at  Wj  =  \/d, 
and  that  as  w  increases  past  wH  to  unity,  5  decreases  mono¬ 
tonically  to  zero;  we  therefore  need  discuss  only  w  <  u\t .  The 
square  root  in  ( A4)  vanishes  if  —  /  =  e/0 ~  1.202;  the  asso- 
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ciated  stationary  point  w„  =  —  (5  +  t)/bt  =  2/ 
(5  —  d„)~ 0.5266  is  a  point  of  inflection  (S'  =  S "  =  0), 
after  which S increases  to  oc  at  wd  =  l/rf„=;0.832.  lfd>dlt, 
then  ( A4)  has  no  real  roots,  and  5  increases  monotonically 
from  0  to  oo  as  w  increases  from  0  to  it',,.  For  1  <  d  <dw  both 


roots  of  (A4)  correspond  to  extrema  (iv¥=wA, 
w_  =u\  )  with  values  satisfying  wA  <ivv  <w,,  =  \/d. 
Thus,if</  =  1.2,  then  there  is  a  local  maximum  at  w  A  =  1/2 
followed  by  a  local  minimum  at  tt'v  =  1/1.8  and  a  singular¬ 
ity  at  iv,,  =  1/1.2.  If  d=  1.1.  then  wA  ~0.38,  wv  -0.81, 


andu’j  —0.91.  If d  decreases  to  1,  then  w.  decreases  to  1/3, 
and  u\  and  u\,  approach  1;  in  the  limit  we  obtain  S,  as 
before. 

The  analogous  two-dimensional  E  of  (64)  in  terms  of 
t  =  c  —  1. 

£  =  u-[2  +  (t-  1 )«’ ]/2(  1  -«•)'.  ( A6) 

corresponds  to  hard  convex  disks  for  t  -  1.  Here,  t  =  1  gives 
E ,  for  circles,  and  t  —  —  l  gives  £,  as  for  ( A  l ) .  The  first  two 
derivatives  with  respect  to  w  are 

E '  —  ( 1  +  tw)/(  1  -  if) 

E"  =  (3  +  1  +  2tw)/(  1  -  if)*.  ( A7) 

However,  although  E  increases  monotonically  to  x  with 
increasing  if  for  t  ,•  —  1 .  the  behavior  for  —  t  =  d>  1  is  quite 
different  than  for  (Al).  The  present  ie(/  =  \/d  =  —  \/t 
(corresponding  to  £'  =  0.  and  E  "  negative)  represents  a  lo¬ 
cal  maximum  £  d  /2(d  —  I ):;  as  if  increases.  £  goes  to 
zero  at  if  —  2/(d  +-  1 ),  and  then  approaches  —  x  asm— 1. 
The  corresponding 

S -- u-(  1  —  if)  V(  1  +  fit’)  (A8) 

is  stationary  at 

if  .  |  -  2  ±  (4  -  3/) 1  -  J/3/ .  (A9) 

for  which 

S"  -  —  (  I  -  if  )’(4  •  6 tw)/(  I  +  nc).  if  =  if  .  . 

(A  10) 

For  all  —  t  —  d >  1,  we  see  that  ,S  is  singular  at  u\,  =  l/d. 
but  that  in  distinction  to  (A3)  the  present  function  changes 
sign  for  if  >  u\,  to  approach  0  from  x  ;  we  restrict  consi¬ 
deration  to  if  <if,;.  The  square  root  in  (A9)  vanishes  at 
-  t  d„  =  4/3.  and  if,,  =  1/2  is  a  point  of  inflection  after 
which  S  increases  to  x  at  ic,;  =  3/4.  If  d >4/3.  then  ( A9) 
has  no  real  roots  and  S  increases  monotonically  from  0  to  x 
as  if  increases  from  0  to  u\,.  There  are  two  roots  for 
1  <«/<  4/3;  e.g  ,  if  1.2,  then  if  .  —0.38,  if  —0.73,  and 
if,,  -0.83.  If  </—  1 .  then  S  — S’,  as  before. 
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